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We theoretically study the stability of the solidified 2nd-layer 3 He at 4/7 of the lst-layer den- 
sity adsorbed on graphite, which exhibits quantum spin liquid. We construct a lattice model for 
the 2nd-layer 3 He by taking account of density fluctuations on the 3rd layer together by employ- 
ing the refined configuration recently found by path integral Monte Carlo simulations. When 
holes are doped into the 4/7 solid, within the mean-field approximation, the density-ordered 
fluid emerges. The evolution of hole pockets offers a unified explanation for the measured dop- 
ing and temperature dependences of specific-heat anomalies. We argue that differentiation in 
momentum space is a key to understanding the physics and accounts for multiscale thermody- 
namic anomalies in the mono- and double-layered 3 He systems beyond the mean-field level. 
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Layered 3 He systems adsorbed on graphite offer unique 
playgrounds as an ideal prototype of strongly correlated 
Fermion systems in purely two dimensions. At the 4/7 
density of the 2nd-layer 3 He relative to the lst-layer den- 
sity, 3 He atoms are solidified to form a triangular lattice, 
which is referred to as the 4/7 phase. 1 The 4/7 phase has 
attracted much attention recently, since specific heat 2 
and susceptibility measurements down to 10 /xK 3 have 
suggested the emergence of a gapless quantum spin liq- 
uid. 

Theoretically, the multiple-spin-exchange model 4 has 
been used for the analysis, 5:6 which, however, has left 
fundamental questions about the 4/7 phase: (1) How to 
realize the gapless ground state? (2) How to explain the 
large saturation field of about 10 T 7 ? To resolve these 
issues, we have pointed out the importance of density 
fluctuations between the 2nd and 3rd layers. 8 By con- 
structing a lattice model to describe the 4/7 solid phase, 
we have shown that strong density fluctuations indeed 
stabilize a gapless quantum spin liquid, which can be 
regarded as that caused essentially by the same mecha- 
nism found in the Hubbard model. 9-12 Furthermore, the 
density fluctuation accounts for the enhanced saturation 
field as observed. 8 

In this letter, we report our analysis of the stability 
of the 4/7 phase on the basis of the lattice model, tak- 
ing account of density fluctuations. The model is derived 
from a refined configuration of the 2nd-layer 3 He recently 
revealed by path-integral Monte Carlo (PIMC) simula- 
tions. 13 Our mean-field (MF) results show that a density- 
ordered fluid emerges when holes are doped into the 4/7 
phase and that the evolution of the hole pockets explains 
measured specific-heat anomalies. Our model gives a uni- 
fied explanation of unusual temperature and doping de- 
pendences of specific heat over a range from the 4/7 solid 
phase to the uniform-fluid phase at low densities. We also 
discuss the validity of the present picture beyond the MF 
approximation and possible relevance to the specific-heat 
anomalies measured recently in the double-layered 3 He 
system. 14 



Our analysis starts from the result of recent PIMC 
simulations, 13 which have revealed a more stable con- 
figuration for the 4/7 phase than that first proposed 
by Elser 1 ' 15 as shown in Fig. 1(a). Here, the open 
circles represent the atoms on the 1st layer, and the 
shaded circles, the locations on the 2nd layer when solid- 
ified. If 3 He ( 4 He) atoms are adsorbed on the 1st layer, 
it forms a triangular lattice with the lattice constant 
a=3.1826 (3.1020) A at the saturation density pi =0.114 
(0.120) atom/A 2 . 1 




Fig. 1. (Color online) (a) Lattice structure of the 4/7 phase of 
3 He shown by recent PIMC simulations. 13 Both lst-layer atoms 
(open circles) and 2nd-layer atoms (shaded circles) form trian- 
gular lattices in the solid phase. The area enclosed by the solid 
line represents the unit cell for the solid of the 2nd-layer atoms 
(see text). The lattice constant of the 1st layer is a. (b) Possible 
stable location of the 2nd-layer atoms are shown by circles on 
top of a a X a parallelogram constructed from the 4 neighbor- 
ing lst-layer atoms, (c) Structure of discretized lattice for the 
2nd-layer model. Lattice points are shown by circles. 



The location of the 2nd-layer atoms is in principle de- 
termined as stable points in continuum space under the 
periodic potential of the 1st layer. In the present treat- 
ment, we simplify the continuum by discretizing it with 
the largest possible number of lattice points kept as can- 
didates of stable points in the solid. To illustrate dis- 
cretization, we cut out from Fig. 1(a) a unit cell of the 
1st layer, namely, a parallelogram whose corners are the 
locations on top of the neighboring 4 atoms on the 1st 
layer as in Fig. 1(b). Possible stable locations of 3 He 
atoms on the 2nd layer are points on top of (1) the cen- 
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ters of the lst-layer triangle (Ai), (2) the midpoint of 
two-neighboring lst-layer atoms (A2), and (3) a nearby 
point of the lst-layer atom (A3). Therefore, we employ 
totally 11 points as discretized lattice points in a paral- 
lelogram, shown as circles in Fig. 1(b). Since a unit cell in 
Fig. 1(a) contains 7 parallelograms, it contains 77 lattice 
points in total, as illustrated as circles in Fig. 1(c). Now 
the 4/7 solid phase is regarded as a regular alignment 
of 4 atoms on 77 available lattice points in the unit cell 
shown in Fig. 1(c). 

For inter-helium interaction, we employ the Lennard- 



(d/r) 12 - {a/rf 



where 



Jones potential Vtj(r) = 4e 

6 = 10.2 K and a = 2.56 A. 16 A more refined Aziz 
potential is expected to give similar results under this 
discretization. The interaction between 3 He atoms on 
the 2nd layer is given by Hy = VijriiTij (n^ is the 
number operator of a Fermion on the «-th site) with Vij 
taken from the spatial dependence of Vl j (r) on the lattice 
points. In the actual 3 He system, the chemical potential 
of the 3rd layer is estimated to be 16 K higher than that 
of the 2nd layer. 17 3 He atoms may fluctuate into the 3rd 
layer over this chemical potential difference and it is sig- 
naled by a peak of the specific heat at T ~ 1 K. 18-20 To 
take account of this density fluctuation, we here mimic 
the allowed occupation on the 3rd layer by introducing 
a simple finite cutoff Vcutoff for Vij within the same form 



of Hamiltonian: When Vlj^j) for 



exceeds 



Kutoff, we take Vij = Vcutoff, otherwise Vij = Vlj(?\j)- 
This allows us to take account of the qualitative but es- 
sential part of the possible occupation on the 3rd layer 
by atoms overcoming Vcutoff- 

We consider the spinless-Fermion model on the lat- 
tice H — TJk + Hy where the kinetic energy is given by 
H K = — J2{ij)(tij c l c j + H.C.). We ignore the effects of 
corrugation potential from the 1st layer for the moment 
and discuss it later. By using the unit-cell index s and 
the site index I in the unit cell, we have = r s + 17. Af- 
ter the Fourier transform, c, = c s ,i = J^k c k,ze lk ' rs /VAu, 
the MF approximation with the diagonal order parame- 
ter (rik 1) leads to 
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where the inter-atom interaction is expressed as Vij — 



In i 



V lm {s') with r s 



r s — r t . Then, we have the 
MF Hamiltonian Hmf = Hk + H^f F . By diagonalizing 
the 77 x 77 Hamiltonian matrix for each k, we obtain the 

energy bands H M f = J2kJ2iLx B K k ) c k,i c k,'- 

For the kinetic energy, several choices of £y are ex- 
amined. As noted in ref. 8, if to is determined so as to 
reproduce the total kinetic energy E^ lMC , the main re- 
sult below, measured in the unit of K, is quite insensitive 
to the choice of Uj . Hence, we here show the results for 
tjj = to for the ij pairs up to the shortest-19th r^, i.e., 
for \r,j\ < 2a. The interaction Vij is taken for < 2a, 
since the longer part is ineffective. 8 The recent PIMC 
simulation estimates the kinetic energy of the 2nd-layer 
3 He in the 4/7 phase as E FIMG = 14 K (17 K) for the 
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Fig. 2. (Color online) (a) Parallelogram formed by 4- 
neighbouring lst-layer atoms, (b) The sites for 2nd-layer atoms 
in the parallelogram under corrugation potentials ABi (solid 
circle), AE2 (shaded circle) and AE3 (open circle), (c) Vdtoff 
dependence of the "charge gap" A c for Ai?i=Ai?2=Ai?3=0 
(filled square) and for A£i = -3.0 K, AE 2 = -1.5 K and 
Ai?3 = K (open triangle) in the 3 He/ 4 He/Gr system, and for 
AEi=AE 2 =AE 3 =0 (filled diamond) and for AE 1 = -3.0 K, 
AE 2 = -1.5 K and AE 3 = K (open circle) in the 3 He/ 3 He/Gr 
system. The bandwidth of the 4th band W4 is illustrated for 
AEi=AE2=A_B 3 =0 (filled circle) in the 3 He/ 4 He/Gr system. 
Shaded lines at V^ uto ff = 22.5 K and 60 K are guides for the 
eyes (see text). 



4 He ( 3 He) 1st layer. 13 Thus, we evaluate to by imposing 
the condition, t a {c\cj + H.C.)/(4AT„) = ^ IMC as 

t = 413 mK (502 mK) for the 4 Hc ( 3 He) 1st layer. 

By solving the MF equations for Hmf, we have the 
solution of the \fl x \fl commensurate structure shown 
in Fig. 1(c) with opening the "charge gap" as shown by 
filled squares (diamonds) in Fig. 2(c) when 4 He ( 3 He) is 
adsorbed on the 1st layer. Here, the "charge gap" is de- 
fined by A c = £^ in (k) - Ef ax (k), where Ef (k) denotes 
the minimum or maximum energy of the Z-th band from 
the lowest. We show here the well-converged results for 
large A c and the small-A c region will be discussed later 
for detailed comparison with experiments. We now dis- 
cuss the effects of corrugation potential. PIMC 13, 15 sug- 
gests that the lst-layer atoms make the corrugation po- 
tential on the 2nd layer even within the discretized lattice 
points. It suggests A#i = -3.0 K and AE 2 = -1.5 K 
relative to A.E3 in the notation in Fig. 2(b). This effect 
merely shifts the A c -Vc U toff line toward a larger Vc U t ff, 
as shown by open triangles (circles) for the 4 He ( 3 He) 1st 
layer in Fig. 2(c). Hence, we show our results below for 
AEi = AE 2 = AE 3 = and for 4 He adsorbed on the 1st 
layer as a representative case. 21 We note that Fig. 2(c) 
indicates that the 4/7 phase is more stable when 4 He is 
adsorbed on the 1st layer rather than the 3 He 1st layer. 

The energy band ^(k) for Vcutoff = 60 K at n = 1 
with n = J2iliY,k n Ki/( 77N u) = 4n /77 and its con- 
tour plot in the folded Brillouin zone 22 are shown in 
Figs. 3(a) and 3(b), respectively. When holes are doped 
into the 4/7 phase, in our MF results, the Fermi surface 
appears at the 4th band as hole pockets for no = 0.99 
(Fig. 3(c)) and n Q = 0.97 (Fig. 3(d)) with the density 
order retained. 

This evolution of hole pockets is reflected in thermody- 
namic quantities: As holes are doped into the 4/7 solid, 
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Fig. 3. (Color) (a) The 4th energy band in the folded Brillouin 
zone for —ir/a < k x < ir/a and —ir/a < k y < ir/a with a = 
2^/1 a 22 for V cuto g = 60 K at no = 1. Contour plot of the 4th 
band for (b) n = 1, (c) n = 0.99 and (d) n = 0.97. Hole 
pockets are represented by white regions in (c) and (d). 



a remarkable peak in the specific heat C(T) develops 
at temperatures lower than the density-order transition 
temperature T c , as shown in Fig. 4(a). At the same time, 
C{T) at temperatures right below T c decreases. The tem- 
perature T* at which C(T) has the highest peak for 
T < T c increases as Uq decreases from 1, as indicated 
by arrows in Fig. 4(a). This tendency was actually ob- 
served in the past 18 and recent 20 measurements: The A- 
like anomaly at T = T c ~ 1 K simultaneously with the 
hump at T ~ 40 mK in C(T) observed for n = 0.96 18 
and for no = 0.971 20 is an indication of the coexisting 
density order and fluid. Since our model H is the spinless- 
Fermion model, the entropy per site S(T) at the high- 
temperature limit is given by A;b1ii2. When T decreases, 
S(T) sharply drops at the density order T c and the re- 
maining entropy is released at T ~ T* < T c as shown in 
Fig. 4(b). 

In experiments, sharp decreases in S(T) at much lower 
temperatures, i.e., T ~ 0.2 and 1 mK, are observed for 
0.9 < no < 1.02, 2,20 which result in a double peak in 

C(T). This spin entropy is ignored in the present model. 
These sharp drops correspond to the release of the spin 
entropy at the energy scale of the spin exchange interac- 
tion J. 23 Under hole doping, the double peak at T ~ 0.2 
and 1 mK is suppressed, whereas the hump C(T*) 20 at 
40 mK grows. This is naturally understood by the sup- 
pression of the density order simultaneously with an in- 
crease in the fluid contributions induced by doping. Actu- 
ally, a double-peak structure in C (T) around T ~ J has 
been shown by our exact-diagonalization calculations 8 
for the density-ordered solid phase in a minimal lattice 
model, which was introduced to mimic the 4/7 solid. 

To further clarify the origin of T*, we calculated the 
density of states for Vcutoff = 60 K at uq — 1, as in 
Fig. 5(a). Since the energy gap opens in the 4/7 phase, 
A c - 1 K, the hump structure in C(T) at T ~ T* < A c 
should arise from the 4th band. As holes are doped into 
the 4/7 phase, the chemical potential shifts to lower ener- 
gies inside the 4th band, as shown by vertical lines in the 
inset of Fig. 5(a). In Fig. 4(a), we plot the energy differ- 
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Fig. 4. (Color) (a) Specific heats for no = 1.0 (light- blue dia- 
mond), 0.99 (green circle), 0.98 (blue open triangle), and 0.97 
(red square) for V cuto ff = 60 K. W 4 op — fi for no = 0.99 (green 
dotted line), 0.98 (blue dashed line), and 0.97 (red line), (b) En- 
tropy per site for no = 1.0 (light blue), 0.99 (green), 0.98 (blue), 
and 0.97 (red) for V cu toff = 60 K. Arrows indicate the tempera- 
tures T* at which C(T) has the highest peak for T <T C . 



and the chemical potential fi, Wl° p — fi, by vertical lines. 
We see that Wl op — (i is located at the central position 
of the hump of C(T) for each n . This confirms that the 
characteristic energy of the fluid in the density-ordered- 
fluid phase expressed by Wl° p — fi corresponds to T*. The 
filling dependences of T* (open square) and W^ op — /i 
(open diamond) are shown in Fig. 5(b). 

To make a detailed comparison with experiments, we 
extrapolate the "charge gap" A c by the least-squares fit 
of the data for 35 K < Vcutoff < 45 K assuming the 
form J2n=o an Kutoff ' 24 wmc h is shown by the dashed 
line in Fig. 2(c). We also extrapolate the width of the 
4th band, W 4 = £^ ax (k) - Ef n (k) by the least-squares 
fit of the data for 35 K < V cut off < 80 K (filled cir- 
cles in Fig. 2(c)) assuming the form J2n=o ^"Kutoff an< ^ 
the result is represented by the dashed line in Fig. 2(c). 
Since the "charge gap" of the 4/7 phase is expected 
to be A c ~ 1 K, 8,18 the corresponding Vcutoff is esti- 
mated to be 22.5 K, which is not inconsistent with the 
chemical potential difference between the 2nd and 3rd 
layers, 16 K 17 (the gray line in Fig. 2(c)). Then, W4 
at Vcutoff = 22.5 K is evaluated to be 174 mK. Since 
W A (Vcutoff = 22.5 K)/W 4 (Vutoff = 60 K) ~ 0.43, the 
density of states of the 4th band at V cu toff = 22.5 K is 
inferred to be enhanced 1/0.43 times more than that at 
Kutoff = 60 K, as shown in the inset of Fig. 5(a). Then, 
lV 4 top - fi and T* at V c utoff = 22.5 K are estimated to be 
0.43 times smaller than those at V cu toff = 60 K, which are 
shown by filled diamonds and filled squares, respectively, 
in Fig. 5(b). The slope of the resultant T* is evaluated to 
be T* ~ 4005 mK with 6 = 1 — no, which is quite consis- 
tent with the experimental observation T* ~ 430<5 mK. 20 
This analysis shows that T* may be regarded as the ef- 
fective bandwidth of the holes by quantum-mechanical 
zero-point motions in the solid, which substantiates the 
hypothesized idea of zero-point vacancy. 2 *- 1,25 

We note that the filling dependence of the entropy 
S{T) for T* < T < T c (for example, S(T = 100 mK), 
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Fig. 5. (Color) (a) Density of states for the 1st (blue), 2nd (pink), 
3rd (green), 4th (red), and 5th (light brown) bands for V cuto fi = 
60 K at no = 1. The inset shows the enlargement of the 4th band 
with W^ 013 (gray line) and the chemical potential u for no = 
0.99 (dotted line), 0.98 (dashed line), and 0.97 (solid line), (b) 
Filling no dependence of T* (square) and W^° p — /i (diamond) for 
^cutoff = 60 K (open symbols with dashed lines) and V cuto ff = 
22.5 K (filled symbols with solid lines). 



not shown) shows nearly the same no dependence as 
S = — nolnno — (1 — no)ln(l — no), as experimentally 
observed. 20 This implies that the entropy for distribut- 
ing NS holes in the iV-site system can be accounted for 
by the hole-pockets contribution in the density-ordered- 
fluid phase. 

We note that a broad shoulder structure of G (T) in the 
uniform fluid phase for n$ = 0.3 (not shown) evolves into 
the C(T*) hump as no increases, which finally shrinks to- 
ward no = f as in Fig. 4(a). A similar G(T) evolution was 
observed in the layered 3 He system on the two 4 He lay- 
ers adsorbed on graphite by Neumann et al. when the 
3 He density n increases and approaches n c = 9.9 nm -2 . 
Since G(T) has common features in both systems, 14,20 it 
is natural to interpret the intervening phase observed for 
rii < n < n c with n\ — 9.2nm in ref. 14 as a density- 
ordered fluid stabilized near the density-ordered solid at 
n = n c . This offers a clear and alternative interpretation 
of the sharp transition or crossover around n = 9.2 nm -2 
in ref. 14. 

The density-ordered fluid whose ground-state energy 
is lower than that of the uniform fluid, £dof("o) < 
Efi^l d (n ), is confirmed at least up to 7% hole doping 
for Kutoff = 60 K. The poor convergence of the MF so- 
lution upon further doping prevents us from determining 
the accurate phase boundary. 

Here, we also stress an alternative possibility, namely, 
the emergence of a uniform fluid phase with small Fermi 
surfaces when holes are doped into the density-ordered 
phase, even when fluctuations beyond the MF theory de- 
stroy the density order in the absence of spin order. The 
Fermi surface is defined by poles of the single-particle 
Green function G(k, oS) at a frequency to = 0. 26 While 
ReG changes its sign across a pole through ±oo, ReG can 
also change the sign across a zero defined by G — 0. In 
the solid phase, only the zero surface exists in the k space 
at u> = while only the poles exist at to = for heavily 
doped uniform fluids. When holes are slightly doped, the 
reconstruction of G(k, oS) yields the interference between 



the zeros and the poles, which creates the resultant Fermi 
surface with the coexistence of zeros and poles. Since the 
interference has a significant k dependence at u> = 0, 
small Fermi pockets may appear after the truncation of 
the large Fermi surface 27 even when the density order is 
destroyed upon hole doping in the absence of the spin or- 
der. It is remarkable that not only doped Mott insulators 
but also doped density-order phases show such differen- 
tiation as observed in the 2D 3 He system. This differenti- 
ation in k space can also be the origin of the peak (small 
cusp) at T - J and the hump (peak) at T ~ T* in G(T) 
for no < l 20 (n < n c 14 ), since the former and latter are 
attributed to the contributions from the truncated and 
remaining parts of the Fermi surface in k space, respec- 
tively. A multiscale measurement of G(T) ranging from 
1 mK to 1 K is desired to resolve this issue in layered 
3 He systems. 
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